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We review recent experimental and theoretical developments in the study of the 
sea quark structure of the proton. In the light quark sector, we analyse the recent 
pp and pD Drell-Yan data from the E866/NuSea experiment at Fermilab, and 
their implication on the d/u asymmetry in the proton. The current status of the 
strange content of the proton, including the possible difference between strange and 
antistrange quark distributions and strangeness form factors, is updated. Finally, 
we point out the implications of the possible non-symmetric charm and anticharm 
distributions in the nucleon for HERA event rates at large x and . 



1 Introduction 

The sea of the proton is extremely fertile ground for the study of non-perturbative 
QCD dynamics and its relation to the substructure of hadrons. In particular, 
sea quark asymmetries almost universally signal the presence of interesting 
non-perturbative phenomena, on the background of a perturbative QCD land- 
scape which in comparison is relatively flat and homogeneous. 



2 Light Antiquark Asymmetry 

The E866/NuSea Collaboration recently measured the spectrum ofji^/i^ Drell- 
Yan pairs produced in pp and pD collisions at the FNAL TevatronEl, which has 
for the first time enabled the shape of the d/u ratio to be mapped out over a 
large range of x. The relatively large asymmetry found implies the presence 
of non-trivial dynamics in the u and d sector of the proton sea which does 
not have a perturbative QCD origin. The novel and unexpected feature of 
the E866 data is that the d/u asymmetry peaks at rather small values of x, 
X ~ 0.15, and drops quite rapidly to unity by x ~ 0.3. 

The simplest and most obvious source of a non-perturbative asymmetry 
in the light quark sea is the chiral structure of QCD. From numerous studies 
in low energy physics, including chiral perturbation theory, pions are known 
to play a crucial role in the structure and dynamics of the nucleon, and there 
is no reason to believe that the long-range tail of the nucleoli should not also 
play a role at higher energies. As pointed out by Thomas □, if the proton's 
wave function contains an explicit Tr^n Fock state component, a deep-inelastic 
probe scattering from the virtual tt"*", which contains a valence d quark, will 
automatically lead to a J excess in the proton. This is the essential physical 
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idea behind these expectations, and has been used to address not only the d/u 
asymmetry (see Rem and peferences therein), but also SU(3) flavoiir symmetry 
breaking in thejproton well as asymmetries in the strangea and heavier 

flavour sectors Bd, as discussed below. 

The essential ingredients of the meson cloud model are the meson-baryon 
distribution functions, fMB{y), which give the probability to find a meson, M, 
in the nucleon carrying a fraction y of the nucleon's light-cone momentum. In 
a hadronic basis the only parameters on which these splitting functions de- 
pend are the hadronic vertex functions, or form factors at the nucleon-meson- 
baryon vertices, which are characterised by an effective momentum scale (cut- 
off) Amb- From previous studies we know that the dominant contribution 
comes from the ttN component of the nucleon, though the ttA configuration 
turns out to play a crucial role here also. In Fig. 1(a) we show the ttN and ttA 
momentum distributions for (dipole cut-offs) Attat — 1 GeV, A^^-a = 1-5 GeV. 
The relative magnitudes of these are taken from a comparison of the axial form 
factors for the nucleon and for the A^-A transition, which strongly favours an 
A'-A axial form factor that is significantly harder than that of the nucleon. 
The harder ttA distribution at large y turns out to be quite important phe- 
nomenologically for understanding the E866 data. 




Figure 1: (a) ttN and ttA momentum distribution functions, with dipole form factor cut- 
offs Attjv = 1 GeV and A^a = 1-5 GeV. (b) Contribution from pion sea to total u + d, 
with A.,riV = ^irA = 1 GeV (largest curve), A^jv = 1 GeV, A^a = Iff GeV (middfc), and 
A^jv = AttA = 1-5 GeV (smallest). The dotted curves are the CTEQ4Q and MRS98Q global 
parameterisations. 

The contributions to the antiquark distribution in the proton can be writ- 
ten as convolutions of the meson distribution functions and the antiquark dis- 
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tribution in the (on-mass-shcU) pion: 



Sqix)= ^ / JiUs{y)'f{x/y), (!) 

where the antiquark distribution in the pion is taken from t:N Dreh-Yan ex- 
periments. The contribution to the total u and d distributions from thcLnon- 
perturbative cloud is shown in Fig. 1(b), compared with the CTEQ4Q and 
MRSQSU parameterisations. While at small x the calculated distributions lie 
safely below the parameterisation, at large x the pion cloud already saturates 
the total sea with cut-offs Ajrjv ^ 1 GeV, A^a ^ 1-5 GeV — although one 
should add a cautionary note that the antiquark distribution at large x is not 
determined very precisely. 

Because the meson cloud model is at most a model of part of the non- 
perturbative sea, it can only be reliably applied to describing the non-singlet 
d — u distribution. To reconstruct the ratio d/u from the calculated difference 
we use, following E866, the total d + u from the CTEQ4 parameterization 
□ as input. With the above form factor cut-offs, one can get a good fit to 
the large- a; data □ stemming from the cancellation of some of the d excess by 
the ttA component. On the other hand, at a; < 0.2 the asymmetry is now 
underestimated somewhat with the hard ttA component □. This suggests that 
there may be room for other mechanisms which could account for the missing 
strength. 

Going beyond explanations involving meson clouds, one can also inves- 
tigate the possibility that the bare nucleon itself could be asyiumetric with 
respect to u and d. As suggested long ago by Field and Feynmanll3, the Pauli 
exclusion principle can contribute to the asymmetry on the basis of the u and 
d valence quarks being unequally represented in the proton, thereby affecting 
the likelihood with which qq pairs can be created in different flavour channels. 

In a simple model in which the nucleon is considered to be composed of 
3 quarks in the ground state, insertion of qq pairs split from the incoming 
virtual photon leads to a ratio of antiquarks in the proton which must satiafc 
d : u = 5 : 4. More quantitative estimates based on the MIT bag model Eil 
showe that the normalisation, A^**"'', of the d — u difference arising from Pauli 
blocking could be as large as 25%. Phenomenologically, one can parameterise 
this contribution as (d — u)^''^'' = A^'*"''(a -I- 1)(1 — x)", where a is some 
large power. Because the E866 data implies a softer asymmetry than typical 
global fits of total sea quark distributions would give, empirically the power a 
uld be > 10 rather than the 5-7 that has been common for the total q fits 
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The Pauli effect will produce an excess of d over u over the whole range 
of a:;, so that it cannot lead to any cancellation of the large-x asymmetry. 
To be consistent with the trend of the large- a; data, especially for the d/u 
ratio, one needs therefore to keep the irNN contribution softer than the ttA^A. 
Taking the ttN and ttA contributions calculated with A^rAr = 1 GeV and 
AttA = 1-5 GeV as above, we show in Fig. 2 the combined effects of pions and 
antisymmetrisation. For the latter the exponent a — 14, and the normalisation 
is A^*^"'' !v 7%, which is at the lower end of the expected scale but consistent 
with the bag model calculations EJ. Together with the integrated asymmetry 
from pions, A^ ^ 0.05, the combined value A = A"" + A^''"'' « 0.12 is in quite 
reasonable agreement with the experimental result, 0.100 ± 0.018 from E866 
and 0.148 ± 0.039 from NMC. 




Figure 2: Contributions from pions with A,rJV = 1 GeV and A^a = 1-5 GeV (dashed) and 
from antisymmetrisation (dotted) to the (a) d — u difference and (b) d/u ratio, and the 
combined effect (solid). 



3 How Strange is the Nucleon? 

There has been a lot of discussion recently about strange matrix elements of 
the nucleon E3. Interest in this subject was largely generated by the deep- 
inelastic scattering experiments with polarised targets at CERN and SLAG, 
which implied a large polarised strange quark distribution in the proton. At 
about the same time a measurement of the elastic neutrino-proton scattering 
cross section at lower values of had also suggested a non-zero value for the 
strange axial vector form factor of the proton. 
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In response to these observations, other processes were sought in which 
traces of strangeness in the nucleon could be detected, such as parity-violating 
filfictron scattering. Some of these have since been performed (at MIT-Bates 
E3), while others (at Jefferson Lab) will soon provide valuable data on the 
strangeness radius and magnetic moment of the nucleon. 

In the DIS regime, the difference between strange and antistrange quark 
distributions in the nucleon has recently come to prominence again-syith the 
availability of new v and i> DIS data from the CCFR collaboration Ej, which 
were analysed for a possible non-zero s-s difference. Perturbative QCD alone 
would be expected to produce identical s and s distributions, while any asym- 
metry would have to be non-perturbative in origin. Such an asymmetry arises 
naturally in a meson cloud picture of the nucleon, where the strangeness of 
the nucleon is assumed to be carried by the kaon-hyperon components of the 
physical nucleon. The s and s quarks therefore have quite different origins in 
this model. Because the s quark originates in the A, its distribution is like 
that of the u quark in the proton, roughly ~ (1 — xf^ at large x. The s in the 
kaon, on the other hand, is much harder, ^ (1 — x). Since the Kh. distribution 
function fxhiy) is fairly symmetric around y = 1/2, upon convoluting the 
momentum distribution with the and distributions, this asymmetry is 
largely preserved, leading to harder s distributions at large x in the proton. 

The asymmetry in the kaon cloud model ±ijrns out to be very small, and 
broadly consistent with the CCFR experimentEJ within the given errors for not 
too large values of A/^-a, Fig. 3(a). To obtain the difference s-s-pye have used 
the absolute values of s -I- s from the parameterisations of RefsLiQ. We should 
point out, however, that there exists some controversy regarding the overall 
normalisation of the deep- inelastic neutrino data from which the strange quark 
distribution was extracted, resulting from an apparent inccpsiatency between 
the neutrino data and data on inclusive charm productionllj'E3. In addition, 
the CCFR data were collected with Fe nuclei targets, so that one needs to 
consider possible nuclear EMC corrections in the data analysisE3 before making 
any definitive conclusions about s and s. 

Within the same formalism one can also calculate the strangeness form 
factors of the nucleon at low . The strangeness (Sachs) radius (defined 
in terms of the strange electric form factor) is found to be very small and 
negative, in the vicinity « —0.004 — > —0.008 fm^ for KNA vertex function 
cut-offs of Aka ~ 0.7-1.3 GeV. The strangeness magnetic form factor, on the 
other hand, suffers from spurious contributions arising from the breaking of 
Lorentz covariance due to the use of the. impulse approximation on the light- 
cone (where the calculation is perfopiiedH) . Once these are removed, according 
to the prescription outlined in Refu-3, the strange magnetic form factor turns 
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• SAMPLE data 



cloud 
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Figure 3: (a) Strange— antistrange quark distribution asymmetry in the nucleon. The soUd 
hues correspond to the asymmetry calculated for Ax a = 0.7 GeV (smallest asymmetry), 1 
GeV and 1.3 Ge\<— ilargest asymmetry), while the shaded region represents the uncertainty 
range of the dataEj. (b) Strange magnetic form factor of the proton G'^j^{Q'^) as a function 
of Q2 The shaded region is the kaon cloud prediction, for a = 8— /lower curve) and 3 
(upper curve) GeV. The data point is from the SAMPLE experiment llj. 

out to be smalLand positive, as in Fig. 3(b), consistent with the trend of the 
SAMPLE datatJ, and largely independent of the details of the KNh. vertex 
function. 

4 Intrinsic Charm 

In 1997 the Hllii and ZEUS0 Collaborations at HERA announced an excess 
of events at large x and in e+p neutral current (NC) and charged current 
(CC) deep-inelastic scattering, which prompted numerous speculations about 
whether one has seen evidence for physics beyond the standard model, such 
as leptoquarks or contact interactions. Since then the apparent signal has 
decreased, although an excess still persists. More interesting from the point of 
view of the physics of non-pcrturbative QCD discussed here, the HERA data 
opened up an avenue through which to investigate the structure of the heavy 
quark sea. 

It was suggested that a large enhancement of the cross sections could 
be achieved by slightly increasing some of the quark distributions, such as the 
u quark, at large x. Unfortunately, as pointed out in Refnd, the size of the 
additional u quark contribution necessary to achieve sufficient enhancement 
would significantly overestimate the SLAG large x data. On the other hand. 




6 



the possibility that the charm quark distribution might be enhanced at large 
X, thereby producing a similar effect, has not been ruled out, and indeed, one 
could more easily imagine that the charm quark distribution could be rather 
hard, owing to the large c quark mass. 

The effect of a non-perturbative, or intrinsic, charm component on the DIS 
crosSj-sections at HERA kinematics was recently investigated by Gunion and 
VogtE3 within a raodel of the 5-quark component of tiie nucleon wave function 
on the light-coneE^, with normalisation fixed to 1%E3. After evolving to the 
HERA kinematics, this intrinsic charm distribution, while considerably harder 
than that generated through pQCD, was still too soft to account for the excess 
HERA eventsE3. As an alternative to this intrinsic charm model, the charmed 
sea was taken in Refsul to arise from the quantum fluctuation of the nucleon to 
a virtual + Ac configuration, along the lines of the tt and K cloud models 
discussed above. A natural prediction of this model is non-symmetric c and c 
distributions. 

Because of the large mass of the c quark, one can approximate the c dis- 
tribution in the D~ meson u and the c distribution in the A+ by: {x) w 
6{x — 1), c^^ (x) « 6{x — 2/3). Quite interestingly, the shape of the resulting 
Sc guar If distributions is quite similar to thaiia the intrinsic charm model of 
RefsESO. However, although the model ofc^cj assumes identical shapes for 
the non-perturbative c and c distributions, the meson cloud gives a significantly 
harder c distribution. 

Evolving the charm and anticharm distributions to an average value of 

= 20000 GeV^, the effects on the CC cross section are seen in Fig. 4(a). 
With a 0.5% (1%) intrinsic charm component the CC cross section increases 
by a factor ^ ^3) for 200 < ill < 250 GeV, which is similar to the excess 
observed by H1II3 in this region. 

While the values-|Of corresponding to the x and values are too low 
for the SLAG data Eil to be sensitive to the charm component of the nucleon 
wave function, data are, however, available from the BGDMS Collaboration 
in the region x > 0.6 which are above charm threshold EZI, and can be used 
to provide limits on the size of the intrinsic charm. The effect of adding 
the intrinsic c distribution to the data on the deuteron structure function is 
illustrated in Fig. 4(b). With the addition of the_D.5% contribution there does 
not seem to be any inconsistency with the dataE3 and, indeed, the agreement 
is slightly improved at x — 0.75. On the other hand, the 1% case may be a 
little too high for comfort. 

Another possibility which could lead to additional enhancement of the 
cross sections at large x would be a larger d quark distribution. The recent 
reanalysis, for example, of the SLAG deuteron data in RefcJ suggested that 
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CC Q^ = 2C000GeV' 
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Figure 4: (a) Ratio of modified to standard DIS model CC cross sections at = 
20000 GeV^, with the modifications arising from 0.5% and 1% additional Sc distributions, 
as well as a modified d quark distribution at large x. (b) Structure function of the deuteron, 
with 0.5% (solid) and 1% (dashed) additions of non-perturbative charm to the global 
fit (dotted) from NMC. 

the valence d/u ratio does not tend to zero as a; ^ 1, as assumed in global 
fits to data 13, but rather Jls consistent with the expectation of perturbative 
QCD, namely d/u — > 1/5E3. The effect of the modified d quark distribution 
is comparable to that of the 0.5% anticharm scenario, as seen in Fig. 4(a). 
The fact that all standard sets of parton distributions assume that d/u ^ as 
X — > 1 means that there is a possible source of systematic error in the modeling 
of "background" events which should be accounted for. 

5 Conclusion 

The high precision deep-inelastic, Drell-Yan and other experiments in the last 
few years have only begun to unravel the rich substructure of the proton sea. 

We have, for the first time, at our disposal important new data from the 
E866 Collaboration which map out the a;-dependencc of the d/u asymmetry in 
the proton. Most importantly, the E866 results confirm the recent observations 
that the d and u content of the proton is not symmetric, while an interesting 
new feature of the data is the relatively fast downturn in the d/u ratio beyond 
X ^ 0.15, which drops rapidly back to unity by x 0.3. The current evidence 
from large x clearly indicates the necessity of a ttA component in the nucleon 
wave function, one which is harder in momentum space than the ttN compo- 
nent. Consistency with data for the sum of d and it at a; > 0.2 requires that 
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both the ttNN and ttNA form factors be relatively soft, making it difficult 
to avoid underestimating the E866 asymmetry at intermediate x, and leaving 
room for other effects, such as the Pauli exclusion principle, to make up the 
difference. Along the lines of previous estimates of the Pauli effect, we find the 
contribution to the d — u difference from antisymmetrisation to be significant 
in magnitude, and particularly important small at x. Our results suggest that 
the best description of the E866 data is that in which chiral symmetry and 
antisymmetrisation plfiK roughly equal roles — consistent with the findings 
of the earlier analysis E3 of the NMC data for F2 — i^J*. We note, however, 
that it would be helpful to have more data at large x, where the error bars 
are largest, to verify the downward trend of d — u, and to further explore the 
possible discrepancy between the Drell-Yan and NMC data. 

Using the same meson cloud framework, one can estimate the asymmetry 
between the s and s quark distributions in the nuclccp. which has been studied 
by the CCFR Collaboration in v charm production til. The magnitude of the 
s-s difference turns out to be very small, consistent with current experimental 
errors, with the s distribution slightly harder than the s. More definitive 
conclusions will be reached with more statistics on the charm production data, 
and the apparent discrepancy between the inclusivPud aep- inelastic muon and 
neutrino data and those on cc production resolved EJlla^t£l. 

Finally, the task of identifying a possible intrinsic charm component of the 
nucleon can be advanced by tagging charm final states in J/ip production at 
HERA to measure the charm structure function, or through a thorough com- 
parison of NC and CC cross sections, for both e+p and e^p collisions, to enable 
one to determine whether the large- a; enhancement of the cross sections is due 
to charm or other flavours. Needless to say, we eagerly await further results 
from the HI and ZEUS Collaborations. Any experimental evidence supporting 
the suggestion that the intrinsic charm could have a strong asymmetry would 
mean a revision of current wisdom, and would undoubtedly lead us to a deeper 
understanding of the structure of hadronic systems. 
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